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Abstract 

A tethered surface model is investigated by using the canonical Monte Carlo sim- 
ulation technique on a torus with an intrinsic curvature. We find that the model 
undergoes a first-order phase transition between the smooth phase and the crumpled 
one. 
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1 Introduction 



Recently, it has been growing that the concern with elastic surface models 
of Helfrich and Polyakov-Kleinert [1,2,3,4,5]. A considerable number of stud- 
ies have been conducted on the phase transition between the smooth phase and 

the crumpled one over the past two decades [6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,2' 

Curvature energies play a crucial role in smoothing the surface. According to 
curvature energies, surface models can be divided into two classes; one with 
an extrinsic curvature and the other with an intrinsic curvature. It is also 
possible that both extrinsic and intrinsic curvatures are included in a model 
Hamiltonian. Intrinsic curvature is the one that is defined only by using the 
metric tensor (first fundamental form) of the surface, and extrinsic curvature 
is defined by using the extrinsic curvature tensor (second fundamental form) 
of the surface [37]. Both of the mean curvature H and the Gaussian curvature 
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K defined by using the extrinsic curvature tensor, are considered as an ex- 
trinsic curvature. However, the Gaussian curvature can also be considered as 
an intrinsic curvature, because of the relation 2K = R, where R is the scalar 
curvature defined only by using the metric tensor. 

Intrinsic curvature models were first studied by Baillie et. al. in [38,39,40,41]. 
The shape of surfaces can be strongly influenced by intrinsic curvatures. Re- 
cently several numerical studies have been made on the phase diagram of the 
model with intrinsic curvature [23,24] . It was reported that the model under- 
goes a first-order phase transition between the smooth phase and the crumpled 
phase on a sphere [23] and on a disk [24] . As a consequence the phase structure 
of the model has been partly clarified: the transition can be seen independent 
of whether the surface is compact or non-compact. 

However, the model is not yet sufficiently understood. Remaining subject to 
be confirmed is whether the phase transition and the topology-change are 
compatible in the surface model on compact surfaces. Here compatible means 
that both of two phenomena lead to the same result without depending on 
which phenomenon firstly occurs. If we define the surface model on compact 
surfaces, it is reasonable to consider the topology-change of surfaces. In fact, 
the partition function for the closed string model of Polyakov includes the 
summation over topology. Moreover, a toroidal vesicle and a genus two vesicle 
with two holes can be observed in biological membranes [4]. 

Therefore, it is interesting, even in the context of the tethered surfaces, to 
study whether the first-order transition can be observed on a torus, as one of 
the higher genus surfaces. The phase transition might be seen only on spherical 
surfaces, if the phase transition and the topology-change are incompatible. 

In this Letter, in order to confirm the compatibility we will show that a first- 
order transition can be seen in the model on a torus with intrinsic curvature 
and that the phase transition is identical to that observed in the model on a 
sphere reported in [23] . 



2 Model 

A triangulated real torus is obtained by modifying a triangulated rectangular 
surface of size Li x L2 as shown in Fig. 1(a). Bending the surface and connect- 
ing the sides of length Li, we have a cylinder of length Li. Then the remaining 
two-sides of the cylinder are connected as in the first-step. Thus, we have a 
real torus as shown in Fig. 1(b), which is topologically identical to the surface 
in Fig. 1(a) under the periodic boundary condition. The real torus mentioned 
above will be used in the simulations. 
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(a) 



(b) 



Fig. 1. (a) Rectangular surface of size Li x L2 = 20 x 10, and (b) the real torus 
obtained by connecting the opposite sides of the surface in (a). 

The real torus is, therefore, characterized by the ratio L1/L2. Two kinds of 
tori are used in the simulations: L\jL2 = 1 and L\jL2=^- Figure 1(b) is the 
torus of size = 200 of the first type L1/L2 = 2, where L\ = 20 and L2 = 10. 
Every vertex has a coordination number o" = 6 on the torus. 

The Gaussian tethering potential Si and the intrinsic curvature 5*3 are defined 
by 

Si = Y.{X,-X,)\ ^3 = -Elog(V2vr), (1) 

where is the sum over bond (ij) connecting the vertices Xj and Xj, and 
Si in 5*3 is the sum of angles of the triangles meeting at the vertex i, and J2i 
is the sum over vertices i. 

The partition function is defined by 



Z{a) = J l[dX,exp[-SiX)], (2) 
1=1 

Six) = Si + aSs, 

where N is the total number of vertices, which is equal to Li x L2 as de- 
scribed previously. The expression S{X) shows that S explicitly depends on 
the variable X. The coefficient a is a modulus of the intrinsic curvature. The 
surfaces are allowed to self- intersect, and the center of each surface is fixed in 
the partition function Z{a) to remove the translational zero-mode. 
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We note, as described in [23], that the intrinsic curvature term S?, = — log((5j/27r) 
comes from the integration measure Hi dXiqf [42] in the partition function for 
the model on a sphere, where qi is the co-ordination number of the vertex i. 
The term S^ = — log(5i/27r) becomes minimum when 5i — 27r for all i, and 
hence it becomes smaller on a smooth torus than on a crumpled one. It is also 
exact that the term S^~Yj^{5i — 27r)^ can be minimized on smooth configura- 
tions of the torus. The reason for using the term = — J2i log(^i/27r), as an 
intrinsic curvature on the torus as in the model on a sphere [23], is that 5*3 is 
closely related to the previously mentioned integration measure. 



3 Monte Carlo technique 

We use two groups of surfaces classified by the ratio L1/L2 in the simulations 
as mentioned above. The first is characterized by L1/L2 — 2 and is composed 

of surfaces of size iV = 1800, iV = 3200, A^ = 5000, and A^ = 9800. The second 
is characterized by L1/L2 = 4 and is composed of surfaces of size = 1762, 
TV = 3600, AT = 6400, and iV = 10000. 

The variables X are updated by using the canonical Monte Carlo technique 

so that X' = X + SX, where the small change SX is made at random in a 
small sphere in R^. The radius 6r of the small sphere is chosen at the start 
of the simulations to maintain the rate of acceptance rx for the X-update as 
0.4 <rx < 0.6. 

The total number of MCS (Monte Carlo sweeps) after the thermalization MCS 
is about 1.5 X 10^ in the smooth phase at the transition point of surfaces of 
> 5000, and about 1.2 x 10^ for the smaller surfaces. Relatively smaller 
number of MCS (0.8 x 10^ ~ 1.5 x 10^) is iterated in the crumpled phase, 
because if the surfaces become once crumpled they hardly return smooth. 
This irreversibility was also seen in the model on a sphere [23] and in the 
model on a disk [24]. 

A random number called Mersenne Twister [43] is used in the MC simulations. 
We use two sequences of random numbers; one for 3-dimensional move of 
vertices X and the other for the Metropolis accept/reject in the update of X. 



4 Results 

Figures 2(a) and 2(b) show Si/N against a obtained on the type L1/L2 = 2 
surfaces and on the type Li/L2 = 4 surfaces, respectively. We find from these 
figures that the expected relation Si/N — 1.5 is satisfied. Scale invariance of 



4 



the partition function predicts that S'i/A^ = 1.5. We expect that this relation 
should not be influenced by whether the phase transition is of first order or 
not. 
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Fig. 2. Si/N vs. a obtained from the tori of types (a) Li/L2 = 2 and (b) Li/L2 = 4. 
The size of surfaces can be reflected in the mean square size X"^ defined by 

^^ = ^E(^^-^)', ^ = ^E^^- (3) 



In fact, it is expected that the surfaces become smooth in the limit a ^ oo 
and crumpled in the limit a— >0. In order to see how large the size of surfaces 
of the type Li/L2 = 2 is, we plot against a in Fig. 3(a). Figure 3(b) shows 
those for surfaces of the type Li/L2 = 4. We find a phase transition separating 
the smooth phase from the crumpled one both in Figs. 3(a) and 3(b). The 
transition seems as a first-order one because X'^ discontinuously changes at 
finite a. 

Figures 4(a) and 4(b) show log-log plots of X^ against N, obtained just below 
and above the transition point etc where X^ discontinuously changes as shown 
in Figs. 3(a) and 3(b). The straight line denoted by smooth (crumpled) is 
obtained by fitting X^ in the smooth (crumpled) phase just above (below) etc 

by 



X^ - TVVi^, 



(4) 



where H is the Hausdorff dimension. Thus we have 

2.10 ±0.03, (smooth); 
■ 12.2 ±0.7, (crumpled). 
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Fig. 3. vs. a obtained on surfaces of the type (a) L1/L2 = 2 and the type (b) 
-^i/-^2 = 4. The surfaces in (a) are of size from A'^ = 1800 to A'^ = 9800, and those in 
(b) are of size from iV = 1764 to iV= 10000. 
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Fig. 4. Log-log plots of vs. N obtained just above and below the transition 
point ttc on surfaces of the types (a) L1/L2 = 4 and (b) L1/L2 = 4. The straight 
lines in (a) and (b) are drawn by fitting the data to Eq. (4). 



H{Li/L2^4) = < 



2.17 ±0.02, (smooth); 
14.0 ±0.6, (crumpled). 



(6) 



The results H = 2.10(3) in (5) and H = 2.17(2) in (6) imply that the surfaces 
are almost smooth in the smooth phase. On the contrary, H = 12.2(7) in (5) 
and H = 14.0(6) in (6) imply that the surfaces are highly crumpled in the 
crumpled phase close to the transition point. 

The phase transition can also be reflected in the convergence speed of the MC 
simulations. In order to see the convergence speed, we plot the autocorrela- 
tion coefficient ^(^2) in Figs. 5(a) and 5(b). The autocorrelation coefficient is 
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Fig. 5. The auto-correlation coefficient ^(5*2) in the smooth and the crumpled phases 
on the tori of the types (a) L1/L2 = 2 and (b) L1/L2 = 4. The auto-correlation 
coefficient A{X'^) obtained on the tori of the types (c) Li/L2 = 2 and (d) Li/L2 = 4. 



defined by 



AfQ \ — Si '52 (Tj) 5^2(^-1+1) , s 

Tj+i = Ti + n X 500, n = 1,2, ■ • ■ , 

where 5*2 is the bending energy defined by S2 = J2i{^ ~ cos 9i). The horizontal 
axes in the figures represent 500 xn (n = 1, 2, • ■ •)-MCS, which is a sampling- 
sweep between the samples S2{Ti) and 5'2(rj+i). The coefficients ^(5*2) in Figs. 
5(a) and 5(b) are obtained on surfaces of L1/L2 = 2 and L1/L2 = 4 respectively. 
Prom these figures, we can find that the convergence speed for 5*2 in the smooth 
phase is 10 times or more larger than that in the crumpled phase. 

Figures 5(c) and 5(d) show the autocorrelation coefficient yl(X^) for the mean 
square size X'^. Contrary to the behavior of A{S2), there is no difference in 
A[X'^) between the crumpled phase and the smooth phase. The reason why 
the convergence speed for in the smooth phase is almost identical to the 
one in the crumpled phase is that the phase space volume (C R^), where the 
vertices take their values, is almost the same both in the smooth phase 
and in the crumpled phase. The shape of surfaces in the smooth phase is 
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considered to be almost fixed, although they extend to relatively large region 
in R^. This behavior of A{X'^) is very similar to the one in a model on a sphere 
with extrinsic curvature [31] but very different from the one in the model on 
a disk with intrinsic curvature [24], where A{X'^) in the smooth phase is quite 
larger than that in the crumpled phase close to the transition point. Therefore, 
we understand from ^(^2) and A{X'^) in Figs. 5(a) - 5(d) that the toroidal 
surface can fluctuate only locally, and there is no large-deformation of the 
surface in the smooth phase. It is quite likely that the shape of surface in the 
smooth phase is almost unchanged in the model on compact surfaces, and it 
largely changes in the model on non-compact surfaces so that the phase space 
volume may be relatively large. 

The phenomenon of critical slowing down can not be seen both in A{S2) and 
in A(X'^) at the transition point. The reason of this is because of irreversibihty 
of the transition. The transition from the crumpled phase to the smooth one 
appears to be irreversible. As mentioned in the chapter 3, if the vertices once 
localize to a small region in at the transition point, they hardly expand to 
be a smooth surface. 
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Fig. 6. S5/N vs. a obtained on the tori of (a) Li/L2 = 4 and (b) Li/L2 = 4. 

The intrinsic curvature energy S^/N is shown in Figs. 6(a) and 6(b), and 
it is clear that S3/N changes discontinuously against a. This indicates the 
existence of a first-order phase transition. The variation of S2/NB against a is 
shown in Figs. 7(a) and 7(b), where A^^ is the total number of bonds, which 
is given by Nb~SN. The discontinuous change of S2/NB indicates that the 
phase transition separates the smooth phase from the crumpled phase. 



Snapshots of surfaces of the type Li/L2 = 4 are shown in Figs. 8(a) and 8(b), 
which are obtained at q;= 1000 and a = 1050. The size of surfaces in Figs. 8(a) 
and 8(b) is of iV= 10000. Sections of the surfaces in Figs. 8(a) and 8(b) are 
shown in Figs. 8(c) and 8(d) respectively. It is easy to see that the surface in 
Fig. 8(a) is crumpled and that the surface in Fig. 8(b) is smooth. Although 
the surface in Fig. 8(b) seems to be a cylinder or a tube, it is almost smooth 



8 




Fig. 8. Snapshots of surfaces of the type L1/L2 = 4 and of size N = 10000 at (a) 
a = 1000, and at (b) a = 1050, and (c) a section of the surface in (a), (d) sections of 
the surface in (b). Figures are shown in the same scale. 



as can be seen in the section shown in Fig. 8(d). We note that the cyhndrical 
surface in Fig. 8(b) is rotationally symmetric only around the longitudinal 
axis of the surface. 
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The appearance of an axis along such cylindrical and smooth surface in Fig. 
8(b) does not imply a spontaneous symmetry breaking, which can be seen in 
a surface model in [30]. In fact, a real torus such as the one shown in Fig. 1(b) 
is also rotationally symmetric only along the axis perpendicular to the hole 
of the torus. In this sense, a symmetry of the model in the smooth phase is 
equivalent to the one of the real torus shown in Fig. 1(b). On the contrary, 
the surface in the crumpled phase is expected to have a higher symmetry just 
like in the crumpled phase of the model on a sphere and on a disk. 



5 Summary and conclusion 

In summary, our attentions were focused on the compatibility of the topol- 
ogy change and the phase transition in a tethered surface model on compact 
surfaces with intrinsic curvature. In order to see this compatibility, we have 
studied the model on two kinds of tori in R^; one is characterized by L1/L2 = 2 
and the other by L1/L2 = 4, where Li and L2 are the length of two sides of 
rectangles which are deformed to be the tori by connecting the sides. Monte 
Carlo simulations were performed on surfaces of size up to A^(= L1XL2) — 1x10"^ 
in both kinds of tori. 

We found that the model undergoes a first-order phase transition separating 
the smooth phase from the crumpled one on both kinds of tori. The smooth 
phase is characterized by Hausdorff dimension if~2, which is identical with 
the topological dimension of surfaces. The surface in the crumpled phase is 
completely collapsed, and characterized by large values of Hausdorff dimen- 
sion. These results indicate that the phase transition on a torus is identical 
with that on a sphere [23] and that on a disk [24]. Thus we can conclude that 
the topology change and the phase transition is compatible in the tethered 
model of surfaces with the intrinsic curvature. Moreover, the autocorrelation 
coefficients A{X'^) obtained on the tori in this paper and on the disk in [24] 
enable us to consider that compact surfaces deform very little in the smooth 
phase, and on the contrary, non-compact surfaces can largely deform in the 
smooth phase. 

An interesting problem to be studied is whether the compatibility is present 
in the model with extrinsic curvature, which exhibits a first-order transition 
on a sphere [23]. It is also interesting to see how the phase transition of the 
intrinsic curvature model is influenced by the fluidity of lateral diffusion, which 
is realized on dynamically triangulated surfaces. 

This work is supported in part by a Grant-in- Aid for Scientific Research, No. 
15560160. 
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